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Summary. We define and study the multiparameter fractional Brownian motion. 
This process is a generalization of both the classical fractional Brownian motion 
and the multiparameter Brownian motion, when the condition of independence is 
relaxed. Relations with the Levy fractional Brownian motion and with the fractional 
Brownian sheet are discussed. Different notions of stationarity of the increments for 
a multiparameter process are studied and applied to the fractional property. Using 
self-similarity we present a characterization for such processes. Finally, behavior of 
the multiparameter fractional Brownian motion along increasing paths is analysed. 

1 Introduction 

The aim of this paper is to give a satisfactory definition of the concept of Mul- 
tiparameter Fractional Brownian Motion (MpfBm) . The definition given here 
is a particular case of the Set-indexed Fractional Brownian Motion studied in 
[HeMe06], but in the multiparameter case, the various stationarity properties 
can be compared. 

In the last decade, two other definitions for the MpfBm appeared in the 
literature (see [He06] for a review of their properties) . Both are problematic as 
extensions of the classical fractional Brownian motion. In this work, we hope to 
persuade the reader that our definition is natural, is the "right" generalization 
of the fractional Brownian motion (fBm) and can be applied directly to real 
applied problems. 
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2 Definition of the MpfBm 

In [HeMe06] , a set-indexed extension of fractional Brownian motion was de- 
fined and some extensions of fractal properties were established. Let A be 
an indexing collection of compact subsets of a metric measure space (met- 
ric d and measure m) satisfying certain assumptions, the Set-indexed frac- 
tional Brownian motion (SifBm) was defined as the centered Gaussian process 
H H = {Bf ; U G A} such that 

V U, V G A; E [Bf Bf ] = X - [m(U) 2H + m(V) 2H - m(U A V) 2H ] , (1) 

where < H < \ and m is a measure defined on the cr-algcbra generated 
by A. 

As the collection A = { [0, t]; t G R+ } U {0} is a particular indexing collec- 
tion, definition (1) provides a multiparameter process which can be seen as a 
multiparameter extension of fractional Brownian motion. We get the following 
definition. 

Definition 1. The Multiparameter Fractional Brownian Motion (MpfBm) is 
defined as the centered Gaussian process B H = {Bf ; t G R+ } such that 

Vs,t€K»;E [Bf Bf ] - \ [m([0, s\f H + m([0, t]f H - m([0, a] A [0, t]f H ] 

(2) 

where m is a measure on R w and He (0, 1/2] is called the index of similarity. 

This definition looks very natural since it relies on a set-indexed process 
and thus, structure of the space R N is only present in indices and not in the 
shape of the covariance function. 

Notice that the definition of the MpfBm depends on the measure m. 

In the particular case of with the Lebesgue measure m, we can explic- 
itly give the covariance between s = (s\, s 2 ) and t = (ii, ti) 

E [Bf Bf ] = \ [( SlS2 r + (ht 2 f H 

- (sis 2 + ht 2 - 2(8! A ti)(a 2 A t 2 )) 2H ] . 

Let us notice that parameter H is restricted to be in (0,1/2], on the 
contrary to standard fractional Brownian motion, in which H is in (0, 1). 

Remark 1. If the measure m is absolutely continuous with respect to the 
Lebesgue measure, the process B ff is almost surely null on the axis. 

Self-similarity is the first property of MpfBm. As a particular case of the 
set-indexed fractional Brownian motion, the multiparameter process inherits 
its properties. It is self-similar of index N.H: for all a G R+, 

{Bf^GR^K^Bf^eR?}. 
where == denotes equality of finite dimensional distributions. 
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3 Comparisons with other multiparameter extensions of 
fBm 

The following two multiparameter extensions of fractional Brownian motions 
are classical. Their definitions rely on a generalization of covariance struc- 
ture of fBm based on euclidian structure of H N . The first definition uses the 
euclidian norm and the second one uses the canonical basis of R N . 

3.1 The Levy fractional Brownian motion 

The Levy fractional Brownian motion (Levy fBm) is defined as the mean-zero 
Gaussian process £$ H = {^^; t <E R+ } such that 

Vs, t e r+ ; e [of a*\ = \ [inr + w H - ii* - «n 

where H G (0, 1). 

The structure of the covariance function of 38 H provides an extension of 
fractional Brownian motion where the absolute value in R + is substituted 
with the euclidian norm of the space . From this point of view, the Levy 
fBm is usually called an isotropic extension of fBm. However, with this sim- 
ple generalization, the process does not seem to be really a multiparameter 
process. 

This simple definition allows to state directly the self-similarity property. 
For all a £ R+, 

{^;i6H^}i{flX;teR?}. 

3.2 The fractional Brownian sheet 

The fractional Brownian sheet is defined as the mean-zero Gaussian process 
B H = {Bf ; t e R+ } such that 

, N 

Vs, t e R+; E [Bf Bf ] = \ JJ + ^ \U 

i=i 

where H= (H 1 ,...,H N ) e (0,l) w . 

In this definition, the euclidian structure of the space H N is strongly 
present in the shape of the covariance function of the fractional Brownian 
sheet. Particularly, this kind of tensor product of standard fractional Brow- 
nian motions along each direction of the canonical basis of R w seems quite 
artificial and lacks generality to be really efficient in concrete applications. 

From the covariance structure of fractional Brownian sheet, the self- 
similarity property can be easily established. For all a e R + , 

{Bf t ;t€R^}®{aS,^]Bf;t€R^} 
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4 Different notions of stationarity 

Stationarity of increments is one of the two characteristic properties of the 
classical fractional Brownian motion. In the framework of multiparameter 
processes, the notion of stationarity can take different forms: 

• Stationarity against translation 

\/h g R+; {X t - X ; t g R+ } = {X t+h - X h - t g R+ } (3) 

• Stationarity in the strong sense 

Vgeg(R N ); {X t -X ; t e R^} ^ {X g(t) - X g(0) ; i e R^ } (4) 
where Q (R") is the set of rigid motions on H N ; see [SaTa94, p. 392]. 

For the next definitions, one needs the notion of the increment of a process 
X on a rectangle D — [s,t], s — (si, . . . , s/v) and t = (ti, . . . , tjv) where s -< t 
(« i <t i ,i = l,...,JV) 

AY(I>)= ]T (-l)^^^^. 

r£{0,l} w 

This definition can be extended to finite unions of rectangles of R+. For 
C = U"=i wnere the Di are rectangles such that 

Vt,j G {i,...,n} Anflj/N^j, 

the increment ZiX(C) is defined by 

n 

AY(C) = ^AX:(A). 

i=l 

This definition is consistent as the previous expression is independent of the 
representation of C . 

• Increment stationarity against translation 

V7i G R+; {AX^; t G R+ } - {% H k]l * e R^} (5) 

• Increment stationarity in the strong sense 

VgeG (R N ) ; {Ay [0 , t] ; t G R+} = {M b(0)ij(()1 ; i £ R? } (6) 

• Measure stationarity (also called Co-incrcment stationarity) 

Vt, Vr >- r' G R^; m([0, r])-m([0, r']) = m([0,t]) => X 4 -X = X T -X T , 

(7) 
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• Increment measure stationarity 

For all finite unions of rectangles C and C, 

m{C) = m(C) AX C = AX C >. (8) 

Notice that, among these 6 properties of stationarity, the first 4 are process 
properties, but the last 2 properties are pointwise properties and depend of 
the chosen measure m. 

The following result summarizes the connections between these different 
stationarity properties: 

Proposition 1. The following implications hold: 

(4) (3) => (5); 
(4) (6) => (5); 
(8) (7). 

From proposition 3.6 and theorem 4.4 in [HcMc06], the following can be 
stated: 

Proposition 2. The MpfBm is Co-increment stationary, but not increment 
measure stationary if H ^ i . 

Let B H be a MpfBm. The increment covariance between two rectangles D 
and D', E[AB H (D) ■ AB H (£>')] can be computed, but the formula is quite 
complicated. 

In the particular case of R+, with the Lebesgue measure and D = D' = 
[s,t], we get: 

E[AB H (D)] 2 = (hh - Sl t 2 ) 2H + (ht 2 - t lS2 ) 2H - (ht 2 - s lS2 ) 2H 

- (sit 2 + hs 2 - 2s lS2 ) 2H + (sit 2 - sis 2 ) 2H + {tis 2 - sis 2 ) 2H . 

We summarize stationarity properties for other definitions of multiparam- 
eter fractional Brownian motion (see [He06], [SaTa94] and [HcMc06] ) . 

Proposition 3. The Levy fractional Brownian motion SS H (H e (0, 1)) sat- 
isfies (3), (4), (5), (6), and the fractional Brownian sheet TB H (H 6 (0, 1) N ) 
satisfies (5). Moreover, if ~B H has constant parameter H in every axis, then 
it satisfies (7). 

5 Characterization 

In Sections 2 and 4, the multiparameter fractional Brownian motion was 
shown to be self-similar and Co-increment stationary. As standard fractional 
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Brownian motion is characterized by its two fractal properties, self-similarity 
and stationarity, it is natural to wonder what are the multiparameter processes 
satisfying the two properties. 

As a particular case of set-indexed fractional Brownian motion, the mul- 
tiparameter fractional Brownian motion satisfies a pseudo-characterization 
property. 

Proposition 4. Let X = {X t ; t G R+ } be a multiparameter process satisfy- 
ing the following two properties: 

1. self- similarity of index a € (0, N/2), 

2. Co-increment stationarity, for Lebesgue measure m. 

Then, the covariance function between s and t such that s -< t is 

E [X s .X t ] = K [m([0, s]) 2a l N + m([0, t]) 2a ' N - m([0, t] \ [0, s]) 2a / N 

Proof. The result simply relies on Proposition 4.1 of [HeMe06], where we 
consider the operation of R + on {[0, t]; t s R-+} such that 

Va>0,VieR^; a.[0, t] = [0, at}. 

In that framework, we have 

Va > 0,Vi e R+; m(a.[0,t}) = a N m([0,t}) 

and then, [i is the function a *—> a N , which is surjective. □ 



A consequence of Proposition 4 is that the fractal properties of self- 
similarity and Co-increments stationarity prescribe the covariance between 
points s and t that are comparable for the partial order -< of R*. Since 
there are non ordered points, we cannot get a complete characterization of 
the MpfBm by the two properties of self-similarity and stationarity. 

A natural question is then, what are the self-similar processes which are 
stationary in the different definitions of Section 4? The following result shows 
that for some choice of stationarity definition, we obtain characterization of 
the Levy fBm. 

Proposition 5. Let H G (0, 1). The Levy fBm is the only Gaussian process 
which is self-similar of index H and stationary in the strong sense (property 

(4))- 

Proof, (cf. [SaTa94, p. 393]) 

It is known (Sections 3 and 4) that the Levy fBm is self-similar and has 
stationary increments in the strong sense. 

Conversely, let X = {X t ; t € R+ } be a Gaussian process such that 
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VaeR+; {X at - teR,^} ( - } {a H X t ; t e R^} 

and 

Vg G Q (H N ) ; {X t - X Q ; t G } ( = } {X fl(t) - X g(0) ; * G R? } 

First of all, considering the canonical basis (ei)i<i<Ar °f R N , and the 
rotation g u that maps ei onto any unit vector w, the slationarity property 
leads to 

E[X u }=E[X gu{ei) -X }=E[X ei ] 
For any s and t in R+, the self-similarity property leads to 

E [X t - X s ] = E [X t _ s ] - E [X ] = \\t - s\\ H E [X ei ] . 
As we also have 

E[X t -X s } = (\\t\\ H -\\s\\ H )E[X €0 ], 

we get E [X t ] = for all (eRf. 

In the same way, we prove that for any s and t in R+ , 

E [(X t X s f] = E [(Xt- S X f] =\\t- sf H E [Xl\ . 

The result follows. □ 

In the fractional Brownian sheet case, several supplementary assumptions 
arc needed to obtain a characterization of the process. Particularly, a null 
value of the process on each axis must be imposed as well as a condition of 
self-similarity for each parameter, when the N — 1 other ones are fixed (see 
[LePo99]). From that point of view, the fractional Brownian sheet has no 
real motivation to be considered, although it satisfies the two properties of 
stationarity and self-similarity. 

6 Projection on flows and regularity 

The notion of flow is the key to reduce the proof of many theorems. It was 
extensively studied in [Iv03] and [IvMeOO]. 

Definition 2. Let S = [a, b] C R. An increasing function f : S — > R^ 
(x < y f(x) -< f(y)) is called a flow. 

The following results, proved in [HeMe06], give a good justification of the 
definition of the MpfBm. 

Proposition 6. Let H H be a MpfBm and f be a flow. Then the process 
(B H y = {B? (4 s, t G [a, b}} is a time changed fractional Brownian motion. 
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However, in general, the projection of a multiparameter process does not 
inherit its different properties. 

Proposition 7. Let f be a flow, and X be a multiparameter process. 

1. If X is a Levy fBm, then (X)f is a classical fractional Brownian motion 
iff f{t) — at where a G R+. 

2. If X is a fractional Brownian sheet, then (Xy is a classical fractional 
Brownian motion iff f is a line parallel to one axis o/R+. 

We conclude this section by giving an interpretation of H parameter. 
Let us recall the definition of the two classical Holder exponents of a 
stochastic process X at to £ R+ : 

• the pointwise Holder exponent 

{\X — X I 
a : limsup sup — < oo > 
p^o s,tets(t , P ) P a J 

• the local Holder exponent 

~ , x f \X t -X s \ } 
Oixito) = sup < a : limsup sup — ; : — < oo > 

{ P^O s ,t€B(t ,p) \t-8\ a J 

Corollary 1. Let B H be a multiparameter fractional Brownian motion with 
self- similarity index H . The pointwise and local Holder exponents of the pro- 
jection (B H )f along any flows f at to £ [0, 1], satisfy almost surely 

<*(B»y(to) = < H otherwise 



~ (j. \ f ®o 



&e(to)-H ifae(t ) < 1 
otherwise 



where is the real function such that 0(t) = m [fit)} (it 6 [0, 1]), and ceg(to) 
(resp. as (to)) is the pointwise (resp. local) Holder exponent of 9 at to- 

Consequently, the H parameter of the MpfBm B ff represents the regular- 
ity of the projection on any regular flow. This fact gives a way to estimate H 
from real data, in the frame of applications. 
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